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Abstract 

Even though the abundance and evolution of clusters have been used to study the cosmological 
parameters including the properties of dark energy owing to their pure dependence on the geometry 
of the Universe and the power spectrum, it is necessary to pay particular attention to the effects 
of dark energy on the analysis. We obtain the explicit dark energy dependent rms linear mass 
fluctuation as which is consistent with the CMB normalization with less than 2 % errors for general 
constant dark energy equation of state, ujq. Thus, we do not have any degeneracy between erg 
and the matter energy density contrast Q^. When we use the correct value of the critical density 
threshold 5 C = 1.58 obtained recently [30, 31] into the cluster number density n calculation in the 
Press-Schechter (PS) formalism, n increases as compared to the one obtained by using S c = 1.69 
by about 60, 80, and 110 % at z = 0, 0.5, and 1, respectively. Thus, PS formalism predicts 
the cluster number consistent with both simulation and observed data at the high mass region. 
We also introduce the improved coefficients of Sheth-Tormen (ST) formalism, which is consistent 
with the recently suggested mass function [36]. We found that changing coq by Awq = — 0.1 from 
luq = —1.0 causes the changing of the comoving numbers of high mass clusters of M = 1O 16 /i -1 M0 
by about 20 and 40 % at z = and 1, respectively. 



1 Introduction 



The formation of the large scale dark matter (DM) potential wells of clusters is solely determined 
by gravitational physics irrelevant to the gas dynamical processes, star formation, and feedback. 
Also clusters are the largest virialized objects in the Universe with their abundance and evolution 
simply related to the linear matter power spectrum. Thus, the abundance of clusters and their 
distribution in redshift should be determined purely by the geometry of the Universe and the 
power spectrum of initial density fluctuations. As a result, the clusters of galaxies provide a 
useful probe of the fundamental cosmological parameters including the investigation into the dark 
energy equation of state ojq, because the linear growth factor D g , the cosmological volume element, 
as well as the primordial scalar amplitude at horizon crossing 5h depend on loq [1, 2]. 

Abundance of rich clusters has been commonly used to constrain the matter power spectrum 
because it is sensitive to the normalization of the power spectrum on cluster scales [3, 4]. The 
amplitude of matter perturbations is sensitive to the presence of dark energy, through the nor- 
malization of the matter power spectrum to the raw large scales probed by the cosmic microwave 
background (CMB) anisotropies [5]. The normalizations of the matter power spectrum obtained 
from both methods can be represented by the rms mass fluctuation on 8/i _1 Mpc scales. How- 
ever, there exists discrepancy in as values resulted from two different normalization methods [6]. 
While erg is almost constant in the former method, it drops rapidly for larger ojq because of the 
increasingly strong integrated Sachs- Wolfe (ISW) effect in the latter. 

There have been numerous papers investigating cluster abundances in the quintessential uni- 
verse [7, 8, 9, 10, 11, 12, 13]. Most of them focus on the influences of dark energy on the background 
evolution and the growth factor. However, we investigate the dark energy effect on the normal- 
ization of the primordial density fluctuation and obtain proper as for the cluster abundance. We 
also use the correct critical threshold density <5 C = 1.58 in our analysis. 

In this paper, we briefly review the effects of quintessence field on the matter power spectrum 
given in Ref. [1]. It has been commonly assumed that there is no significant effects of quintessence 
field on the normalization as because it clusters gravitationally on large length scales but remains 
smooth like the cosmological constant on small length scales due to the relativistic dispersion of 
its fluctuations. However, we directly show the effects of the quintessence field on the as(M) due 
to the change in Sh in the next section. We also show the effects of quintessence field on both 
the background evolution and the growth factor. In Sec. 3, we investigate the resulting changes 
in the cluster abundance. We conclude in Sec. 4. We also consider ujq < — 1 cases, because the 
obtained formalisms can be extended to the general constant ojq models. 
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2 Cosmological Consequences of Quintessence Models 



We briefly review and find the influences of the quintessence field on the cosmological quantities in 
this section. We consider the cold dark matter with the quintessence field (QCDM) in a spatially 
flat universe. We limit our analysis on the constant equation of state of the dark energy ujq which 
is proper for the late time behaviors of quintessence models [14, 15]. We are able to apply these 
solutions to the time- varying ujq by interpolating between models with constant ujq [16, 17]. 



2.1 The Power Spectrum 



The linear perturbation equation for the quintessence field Q = Qo + 5Q is given by 

5Q + 3H5Q + (k 2 + V QQ )5Q 



(2.1) 



where dots mean the derivatives with respect to the cosmic time, V QQ = d 2 V/dQ 2 \ q=q , and h 
is the trace of the spatial metric perturbation [18, 19]. Thus, the Compton wavelength of the 
quintessence field above which it clusters gravitationally but remains smooth on smaller scales is 
determined from the wavenumber kq = 2^/Vqq, where Vqq is given by 
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(2.2) 



where ujq is the equation of state of Q field, fi m (a) and ^de(«) are the energy density contrasts 
of matter and the quintessence, respectively. The linear perturbation of the quintessence field 8Q 
grows only on large scales (k <C fc<g), so the quintessence field clusters and affects the evolution 
of the matter density perturbation 8 m . The Compton wavenumber is determined by two terms 
in Eq. (2.2). The second term is due to the time variation of ujq and disappears when ujq is a 
constant. We investigate the contribution of the second term to fcq when ujq varies slowly. For 
this purpose, we adopt Chevallier-Polarski-Linder (CPL) parametrization ujq = ujq + oj a (l — a) 
[20, 21]. In Fig. 1, the dot-dashed, solid, and dashed lines correspond to (wo,w a ) = (—0.9,0.1), 
(—0.98,0.01), and (—0.8,0.3), respectively. It is less than 30 % for the slowly varying loq cases. 
Thus, we can safely consider only the first term of the wavenumber for those cases. The effects 
of the quintessence field on the matter power spectrum and the time evolution of gravitational 
clustering is parameterized by the shape parameter Tq = kqjh in Ref. [1] and we adopt this 
shape parameter. 
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Figure 1: The ratio of the second term to the first term of of Uq for CPL parametrization values 
(uo,u a ) = (—0.9,0.1), (—0.98,0.01), and (—0.8,0.3) (for dot-dashed, solid, and dashed lines, 
respectively) . 



The linear power spectrum for S m in QCDM models is given by 



P(k,a)=A Q k n *T%(k)\ 



EM 

D(a ) 



(2.3) 



where Aq is a normalization, n s is the spectral index of the primordial adiabatic density pertur- 
bations, Tq(K) is the transfer function, D is the linear growth factor, and the present scale factor 
normalized as ao = 1. Even though the transfer function Tq(/c) does depend on ojq, its change 
happens only on large scale k < 0.01 hMpc -1 . Thus, this correction in the transfer function does 
not affect the value of the rms mass fluctuation. The main effect of the quintessence field is on 
the normalization Aq which can be written as Aq = 27r 2 e)|^(c/-£fo) ns+3 ! where [1] 

2.05 x 10" 5 Q 1 (O m ) Cl+C2lnnm exp[c 3 (n s - 1) +c 4 (n s - l) 2 ] ,with (2.4) 



Sh 
ci 
a 



-0.789|wq 



0.0754-0.211 In \wq\ 



, c 2 = -0.118 - 0.0727w Q , c 3 = -1.037 , c 4 = -0.138 , 



(2.5) 



-ojq) s with s = (0.012 - 0.036o;q - 0.017o; Q i ) ( 1 - n m (a) i 
+(0.098 + 0.029w Q - 0.0850^) lnO m (a) . 

Note that the notation ao = ol{clq) and 5h is the amplitude at horizon crossing. We use 5h = 
2.05 x 10~ 5 to be consistent with WMAP 7 year data, and the cosmological parameters = 0.272, 
h = 0.7, n s = 0.963, Og = 0.0456, and cr 8 = 0.809 from WMAP + BAO + H measurements 
given in Ref. [22]. The rms linear mass fluctuation at top- hat smoothing scale R is given by 

2\ 



8M 



1 



k 2 P(k,a) W(kR) 



dk, 



(2.6) 



M(R, a) 

where the filtering radius R is the Lagrangian radius of a halo of mass M at the present time, 
47rpl4a ) ) 1/3 ' and W(kR) = jj^s(am[kR] - (kR)cos[kR]\ is the top-hat window function. 



R 



3 



0.001 0.01 0.1 1 -11 -10 -0.9 -0.8 -0.7 -0.6 -0.5 

k [h Mpc" 1 ] °> Q 

Figure 2: a) The linear matter power spectrum for different values of ujq = —1.1, —1.0, and —0.8 
(from top to bottom), b) ag values from the COBE normalization and the Wang and Steinhardt's 
approximation (from bottom to top). 

Traditionally, the cluster abundance is used to put a constraint on the dispersion of the density 
contrast at the scale 8 h Mpc, denoted ag. We also denote the mass inside a sphere of radius 
Rg = 8 h Mpc as Mg = 5.95 x 10 14 J1^ 1 /i _1 Mq where M means the solar mass and we use the 
present critical energy density p® Tit = 2.775 x lO 11 M0/i 2 Mpc~ 3 . Now we show the present matter 
power spectra and ag values for different values of ujq in Fig. 2. In the left panel of Fig. 2, 
the dot-dashed, solid, and dashed lines correspond to the matter power spectrum for cjq = —1.1, 

— 1.0, and —0.8, respectively. The differences in the power spectra of ujq = —1.1 and —0.8 from 

— 1.0 are about 7 and 17 %, respectively. The amplitude of the power spectra are comparable with 
the one in Ref. [23]. In the right panel of Fig. 2, we also compare the COBE normalized ag (solid 
line) with the one from Ref. [4] (dashed one). There exists the discrepancy in the dependence 
of ag on ujq between the two different normalizations other than the magnitudes. It is easy to 
understand this as shown in Ref. [6]. For the fixed value of f^, the dark energy dominates the 
cosmic expansion earlier as ujq increases and thus enhancing the dynamics of the gravitational 
potential which results in an increasing integrated Sachs- Wolfe (ISW) effect on large scales. This 
effect is properly shown in the COBE normalized ag only. Also ag obtained from the X-ray cluster 
population as in Ref. [4] suffers from the uncertainties due to the uncertainties in modelling. Thus, 
we use the COBE normalized ag in our cluster abundance calculation. 



Even though the COBE normalized ag is the proper and the accurate one, we need to express it 
as a function of M in the cluster abundance calculation. When ujq = —1, ag is well known [24] 
and we can parameterize cjq dependence of ag by ag(ujQ) = (— wq) - 72+0 ' 36 ''''Qct8(u;q = — 1) with 
less than 2 % errors for —1.1 < ujq < —0.6. For the cosmological parameters adopted from the 
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Figure 3: cr(M, z = 0) for different values of ujq = —1.1, —1.0, and —0.8 (from top to bottom). 



WMAP, we can express 

a(M,z) ~ (- WQ )°- 72 + - 36 ^Q (3.90 - 0.215 log 



,1/ 



A,(*) 



(2.7) 



We show cr(M, z = 0) for different values of ujq in Fig. 3. The dot-dashed, solid, and dashed lines 
correspond to loq = —1.1, —1.0, and —0.8, respectively. 



2.2 Volume Element and Growth Factor 



Friedmann's equation for the QCDM is written as 



H 2 {z) = Hl 



^(l + z) 3 + ^ c (l + z) 3 ( 1+ ^ ) 



Thus, the cosmic volume per unit redshift is also given by 

cdt 



V(z) 



4ird 2 A (z'" 



dz 



{z')dz' , 



(2.1) 



(2.2) 



where d A = Jq hIz 1 ) * s ^ ne an S u l ar diameter distance between redshifts and z. V(z) is the 
proper volume of a sphere of radius z around the observer. We show the dependence of V(z) on 
wq in the left panel of Fig. 4. As wq increases, V(z) decreases. 

The sub-horizon-scale linear perturbation equation with respect to the scale factor a is well 
known [25], given by 

d 2 5 | fdlnH | 3\dS 47rG /0m j = Q 
da 2 \ da a J da (aH) 2 
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z z 

Figure 4: a) V(z) for different values of ujq = —1.1, —1.0, and —0.8 (from top to bottom), b) 
D g {z) for the same values of ujq as in the left panel. 



The exact analytic growing mode solution D g of 5 for any value of the constant ujq is well known 
[16, 26, 27]: 
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(2.5) 



where y = a i fi. In the right panel of Fig. 4, we show the behaviors of the growth factor D g for 
different dark energy models (i.e. for different values of ujq) when = 0.272. The dot-dashed, 
solid, and dashed lines correspond to ujq = —1.1,-1.0, and —0.8, respectively. As ujq decreases, 
D g maintains the linear growth factor proportional to a for a longer time. 
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Figure 5: a) The comoving number density of clusters n of mass greater than M for different 
values of z = 0, 0.5, 1.0, and 2.0 (from top to bottom) when uq = —1 and 5 C = 1.58. The circular 
(z ~ 0) and triangular (0.18 < z < 0.85) dots represent the data from Ref. [33]. b) Errors of n 
when we use the correct threshold density contrast S c = 1.58 instead of 1.69 for different values 
of z = 0, 0.5, and 1.0 (from bottom to top). 



3 Mass Function and Cluster Number 



The mass function f{M) from the Press-Schechter (PS) formalism is related to the comoving 
number density dn of objects in the range dM as dn = (p^ l /M)\dlna/dM\f(M)dM [28, 29]. PS 
theory which relates the comoving number density of virialized objects to their mass is given by 



dn{M,z) = \-^ 



2 o° 
ttM 2 



din a 



dluM 



— exp 



2a 2 



dM, 



(3.1) 



where the critical density threshold <5 C = ^ near is predicted for a spherical overdensity of radius 
R = (3M/47TP™) 1 / 3 and mass M according to the linear theory. However, the PS mass function 
is known to predict too many low mass clusters and too few high mass clusters, as well as too few 
clusters at high z [9]. One remark is that the correct value of 5 C = 1.58 was recently obtained 
independent of the value of wq instead of the well-known value of 1.69 [30, 31, 32]. If we use 
this correct value of 8 C , then the PS formalism shows the improved predictions for the number 
densities of both low and high mass clusters. We show this in Fig. 5. In the left panel of Fig. 5, 
we show the comoving number density of clusters n of mass great than M for ujq = — 1. The solid, 
dashed, dot-dashed, and long dashed lines correspond to z = 0, 0.5, 1.0, and 2.0, respectively. 
The circular (z ~ 0) and triangular (0.18 < z < 0.85) dots represent the data from Ref. [33]. We 
also show the changes of n at high and low masses at different redshifts in the right panel of Fig. 
5. At the low mass M = 10 9 [/i~ 1 M Q ], there is about 2 % decrease in n at present. At the high 
mass M = 10 16 [/i _1 Mq], n increases about 58 % today. Also at high z, n increases about 75 and 
114 % at z = 0.5 and 1.0, respectively. 
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Figure 6: n(> M) from both PS formalism with correct value of 5 C 
formalism with 6 C = 1.69 (solid line). 



1.58 (dashed line) and ST 



Using the correct value of 6 C , we can predict more massive clusters and more clusters at high z. 
However, we may still have too many low mass clusters even by using correct value of 5 C from the 
PS formalism. This is shown is Fig. 6. If we strictly limit the mass of clusters as bigger than 
1O 14 M0, then PS formalism might be good enough for the cluster abundance calculation. Also 
there might be other mechanisms than only gravity for the low mass cluster formations. In any 
case, PS formalism shows the deviation from the simulation at low mass region. Thus, we need 
to consider another popular numerical fit for the differential mass function given by Sheth and 
Tormen (ST) [34, 35], 



/st(o-) 




b5l 
2a 2 




6_c 
a 



(3.2) 



where A = 0.3222, b = 0.75, and p = 0.3 are three parameters tuned to fit with numerical 
simulations. A is fixed by the normalization that all dark matter particles reside in halos. However, 
ST mass function deviates from the simulation results by as much as 40 % at the high mass end 
[36]. Thus, a new fitting function for f(cr) is given in Ref. [36] by adding one extra parameter 
into ST, 



/mod % 

0.333a 011 , b = 0.788o a01 



A\ — exp 

7T 



2a 2 




'drVb 



a 



(3.3) 



where A = 0.333a , b = 0.788a , p = 0.807, and q = 1.795. Recently, there is another mass 
function /Manera which is similar to the original ST formalism [37]. 



However, the ST mass function is known to deviate from the simulation at the high mass end [36]. 
This problem can be cured when we use the correct 5 C . Thus, we introduce a simple but quite 
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Figure 7: The different mass functions as a function of a: /gx, /Manera> /ln 5 an d /mod (from 
bottom to top). 



similar to the results in Ref. [36] by using the original ST formalism: 

0.675 f ' 



/ln(c,2) 



0.32,/?^ exp 



7T 



2a 2 



0.32 



1 + 



0.67(5? 



a 



(3.4) 



The comparison between different mass functions f(cr) is given in Fig. 7. The dashed, long dashed, 
solid, and dot-dashed lines correspond to /st, /ivianera; fhN, and / mo d ; respectively. As we can see 
in this figure, we can produce enough high mass clusters from the simple ST formalism when we 
use the correct 5 r = 1.58. 



We show the comoving number density of the clusters for loq = — 1.0 at different z in the left 
panel of Fig. 8. The solid, dashed, dot-dashed, and large dashed lines correspond to z = 0, 0.5, 
1.0, and 2.0, respectively. In the right panel of Fig. 8, we show the relative errors of the comoving 
number densities for two different models at two different redshifts. For example, the solid line 
represents the relative errors of n between luq = —1.1 and —1.0 at z = 0, |n(> M) UCi= -x,i,z=o~ n (> 
M) WQ= _i.o )2 =o|/ n (> A/)cj Q =-i. 0,2=0 x 100 (%). At low mass M = 1O 9 /i _1 M , the differences of 
n between two different models (cjq = —1.1 and —1.0) are only 0.5 and 1 % at z = and 1, 
respectively. At this mass, the differences of n between ojq = —0.8 and —1.0 are also only 1 
and 4 % at z = and 1, respectively. However, the differences of n between cjq = —1.1 and 
— 1.0 at high mass M = 1O 16 /i _1 M0 are 28 and 48 % at z = and 1, respectively. For the two 
models loq = —0.8 and —1.0, the differences become 58 and 80 % at two different z = and 1, 
respectively. 
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Figure 8: a) The comoving number density for the mass function given in Eq. (3.4) when coq = — 1 
at z = 0, 0.5, 1, and 2 (from top to bottom), b) The relative errors of comoving number density for 
the two different dark energy models between wq = — 1.1 and —1.0 at two different redshifts z = 
(solid) and z = 1 (dashed) and for models between uq = —0.8 and —1.0 at z = (dot-dashed) 
and 1 (long dashed). 

4 Conclusions 

We obtain the dark energy dependence of the mass fluctuation which is consistent with the CMB 
normalization. Even though we limit our analysis for the constant ujq, we can extend our analysis 
for the slowly time varying ujq or use the interpolation of the constant uq to study the general 
time varying ujq. 

We use the correct critical density threshold contrast 5 C to calculate the comoving number 
density of clusters for both PS and ST formalism. We show that PS formalism with this correct 
value of S c can predict the consistent cluster number with both the simulation and the observed 
data at the high mass region. Even though, the PS formalism with this correct value of 5 C still 
predicts too many low mass clusters, this might be due to other mechanisms for the cluster 
formation at the low mass region in addition to gravity. Thus, PS formalism might be good 
enough to explain the cluster abundance. 

We obtain that the dark energy dependence of the comoving number densities for the high 
mass M = 1O 16 /i _1 M0 is about 25 and 40 % at two different redshifts z = and 1. Thus, 
observation will need to be as accurate as this level to probe the property of dark energy from the 
cluster physics. 
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